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Dark matter annihilation is one of the leading explanations for the recently ob- 
served excesses in cosmic rays by PAMELA, ATIC, FERMI-LAT and HESS. Any 
dark matter annihilation model proposed to explain these data must also explain the 
fact that PAMELA data show excesses only in spectrum but not in anti-proton. 
It is interesting to ask whether the annihilation mode into anti-proton is completely 
disallowed or only suppressed at low energies. Most models proposed have negligible 
anti-protons in all energy ranges. We show that the leptocentric U{l)B~3Li dark 
matter model can explain the excesses with suppressed anti-proton mode at low 
energies, but at higher energies there are sizable anti-proton excesses. Near future 
data from PAMELA and AMS can provide crucial test for this type of models. Cos- 
mic 7 ray data can further rule out some of the models. We also show that this 
model has interesting cosmic neutrino signatures. 

PACS numbers: 98.80.Cq ll.15.Tk 11.25.Hf 14.80.-j 



I. INTRODUCTION 

Recently several experiments have reported excesses in cosmic ray (CR) energy spec- 
trum. Last year the PAMELA collaboration reported e"*" excess in the CR energy spectrum 
from 10 to 100 GeV, but observed no anti-proton excess compared with predictions 

from CR physics. These results are compatible with the previous HEAT and AMSOl exper- 



iments (e.g., Ref. 4|) but with higher precision. Shortly after the ATIC and PPB-BETS 
balloon experiments have reported excesses in the e+ + e~ spectrum between 300 and 800 
GeV js, 6]. The ATIC data show a sharp falling in the energy spectrum around 600 GeV. 
Newly published result from FERMI-LAT collaboration also shows excesses in the + e~ 
energy spectrum above the background [3]. However, the spectrum is softer than that from 
ATIC. In addition, the HESS collaboration has inferred a flat but statistically limited + e~ 
spectrum between 340 GeV and 1 TeV which falls steeply above 1 TeV . These obser- 
vational data have generated a lot of excitement as the excesses may be explained by dark 



matter (DM) annihilation or decay with appropriate properties 
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If DM annihilation is responsible for the observed data, a lot of information can be 
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. The mass of the annihilating DM serves 



extracted about DM properties [id, 
as the cut off scale of the spectrum, the lepton spectra must have a cut off energy at 
the DM mass m/). The FERMI-LAT and HESS data would require that the DM mass to 
be around a few TeV. The DM belongs to the weakly interacting massive particle (WIMP) 
category. The fact that PAMELA did not see anti-proton excess indicates that the DM is 
hadrophobic or leptophilic. At least hadronic annihilation modes are suppressed. Also since 
the same DM annihilation rate producing the excesses is related to the annihilation rate 
producing the cosmological relic DM density, 23% of the energy budget of our universe, the 
annihilation rate is constrained. The latter requires that the thermally averaged annihilation 
rate {av) to be 3 x 10~^^ cm^s"^. If this annihilation rate is used to calculate the spectrum, 
it is too small by a factor of 100 to 1000. There is the need to boost up the spectrum with a 
boost factor B of 100 to 1000. Any DM model proposed to explain the data should be able 
to produce the large boost factor too. 

To su ppre ss the cosmic anti-protons leptophilic DM models are proposed by several au- 



thors 
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13| . In this type of models the DM only interacts with leptons in the standard 



model (SM). However, these DM models are hard to be tested since DM does not interact 
with hadrons at all. Another widely considered U{1) extension of the SM is the gauged 
B-L. However, if DM interacts with the SM with the gauged U{1)b-l interaction too much 
anti-protons will be produced and therefore in conflict with PAMELA data. In this work, 
we propose a particle physics DM model, the gauged [15|] U{l)B-zLi DM model, to explain 
the excess data. Here B is the baryon number and Lj is one of the e, and r lepton 
numbers. 
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In this model, the DM is a Dirac fermion with nontrivial U{1) BSLi charge 15j. The gauge 



boson of this U{1) group, Z', mediates DM annihilation into SM particles. This model is not 
a pure leptophhc DM model, like the f/(l)L-_i^ model studied in Ref. [3]. But has larger 
couplings to leptons than quarks. We refer it as the leptocentric dark model. The Z' has 
suppressed couplings to quarks leading to suppressed anti-proton up to the PAMELA energy 
reach, but may lead to excesses at higher energies. With the Z' mass close to two times of 
the DM mass, a large boost factor can be produced through the Breit-Wigner enhancement 
mechanism [ll|, [igI. The shape of the spectrum helps to further constrain the choice of 
the lepton flavor Li of the U{l)B-3Li- 

All DM annihilation models proposed to explain contribute at certain level to 7-rays 
by final state radiation. One must check if the model is compatible with available data. We 
find that for the Lj = r model, the hadronic r decays lead to too much 7-rays in conflict 
with data and only the Li = Le^/j, models can pass the 7 ray data constraint if one takes the 

excess spectrum to fit the ATIC or FERMI-LAT data. It is interesting to note that the Z' 
has the same couplings for charged lepton 1^. and neutrino ul-. Therefore after fitting data 
on excesses in CR, the cosmic neutrinos by DM annihilation from the Galactic center 
(GC) are predicted. The neutrino signals can be tested by experiments such as IceCube. 

In the next section we will first describe how to treat the propagation of charged cosmic 
rays in the Milky Way. Then we present the leptocentric U{l)B~3Li DM model in Sec. III. 
The cosmic e^, p, 7 and neutrino spectra predicted are given in Sec. IV, V, VI and VII 
respectively. Finally we give our conclusions in Sec. VIII. 

II. THE PROPAGATION MODELS 

Before detailed studies of the CRs from the leptocentric U{l)B~3Li DM model, let us 
first briefiy describe the CR propagation model which predicts the positron and anti-proton 
spectra to compare with data, and also 7 and neutrino CRs. 

es propagate diffusively in the Galaxy due to the scattering with 



The charged partic 

random magnetic field 17|]. The interactions with interstellar media (ISM) and interstellar 
radiation field (ISRF) will lead to energy losses of CRs. For heavy nuclei and unstable nuclei 
there are fragmentation processes by collisions with ISM and radioactive decays respectively. 
In addition, the overall convection driven by the Galactic wind and reacceleration due to 



4 



the interstellar shock will also affect the propagation processes of CRs. The propagation 
equation can be written asfisl 



dip 
'dt 



dp 



d_l_ 



d_ 

dp 



p^_|(V.Ve^) 



where ip is the density of cosmic ray particles per unit momentum interval, D^^ is the 
spatial diffusion coefficient, Vc is the convection velocity, Dpp is the diffusion coefficient in 
momentum space used to describe the reacceleration process, p = dp/dt is the momentum 
loss rate, Tf and Tr are time scales for fragmentation and radioactive decay respectively. 
Here Q(x, p) is the source term. For the astrophysical sources they are located within the 
Galactic disk [18|, 119( while for DM annihilation the source is within the whole dark matter 
halo. The source term of DM annihilation is given as 



QiX,E) 



(av) dN 2 



2ml 



(2) 



with (av) the annihilation cross section and the mass of dark matter particle, dN/dE 
the spectrum per annihilation and p(X) the dark matter density at X. 

The spatial diffusion is regarded as isotropic and described using a rigidity dependent 
function D^-j. = PDq j with p the rigidity of the particles. The reacceleration is described 
by the diffusion in momentum space with the diffusion coefficient Dpp, which relates with 
the spatial diffusion coefficient D^x as shown in 20|]. The convection velocit y is cylindrically 
symmetric and increases linearly with the height z from the Galactic plane 18| . Finally the 
energy losses and fragmentations are calculated according to the interactions between CRs 
and the interstellar gas or the interstellar radiation field. 

A numerical method to solve Eg. ([1]) has been developed by Strong and Moskalenko, 
known as the GALPROP model 18|. In GALPROP, the realistic astrophysical inputs 
such as the interstellar medium (ISM) and interstellar radiation field (ISRF) are adopted 
to calculate the fragmentations and energy losses of CRs. The parameters are tuned to 
reproduce all the observational CR spectra at Earth. The GALPROP model can give very 
good descriptions of all kinds of CRs, including the secondaries such as e+, p and diffuse 



7-rays[18, 



21 
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In this work we employ GALPROP to calculate the propagation of CRs. We will adopt 
the diffusion + convection (DC) scenario with parameters given in Ref. 191]. The DC model 
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predicts the e~ , e^, p and p spectra with very good agreement with the data 

The main uncertainty of the propagation model comes from the height of the Galactic 
diffusion region. If the height is changed the diffusion coefficient can be adjusted accordingly 
so that the confinement time of the cosmic particles within the Galaxy and the secondary 
products are kept intact. However, such adjustments affect the contribution from DM 



annihilation greatly 



24| . The reason is that the sources of CRs are located within the 



thin Galactic disk, which is much smaller than the diffusion region no matter how to adjust 
the height. But the contribution from DM annihilation comes from the whole dark matter 
halo. Only when the annihilation takes place within the diffusion region can it contribute 
to the observed cosmic ray fluxes. The particles from the region out of the diffusion region 
escape freely. Therefore the height affects the contribution of DM annihilation greatly. 

The effects by changing the diffusion height are also different for electrons and nuclei 
from DM annihilation. For electrons and positrons the main effects that affect their spectra 
are the energy loss processes, such as the synchrotron radiation or inverse Compton scatter- 
ing with the ISRF. Since electrons/positrons loss energy rapidly, the high energy electrons 



can not propagate for long distance, usually smaller than Ikpc [25|. On the contrary the 
nuclei are mainly affected by the propagation processes and propagate much longer than 
the electrons and positrons. Therefore the diffusion height affects anti-proton contribution 
by DM annihilation much more than positrons. In the following we will try to adjust the 
diffusion height to suppress the anti-proton contribution from DM annihilation while the 
electrons/positrons can fit the PAMELA and Fermi data at the same time. 

In the conventional cosmic ray model the diffusion height Zh is taken as 4kpc, which can 
give best fit to all CR data. Considering the data from diffuse 7-rays the height should be at 
least of about ~ 2kpc. In our calculation we will take another set of propagation parameters 
with the height Zh = 2kpc. All other parameters are adjusted correspondingly to fit the CR 
data. 

For cosmic 7- or neutrino-rays, since they can propagate freely in the Galaxy, one can 
obtain their fluxes by simply integrating the source terms along the line of sight (LOS) at 
any direction. The flux (j){E^^i,) for 7 and neutrino are given by 

1 ^ {av) dN 

^-^l^f J LOS 



where (j){E-^^^, ip) are the 7- and neutrino-ray spectra at the direction ip and dN/dE^^^, are the 
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7" and neutrino-ray spectra per annihilation, respectively. For the emission from a diffuse 
region with solid angle Af2, we define the average flux as 

(i)An = ^J^jWdn. (4) 

In the rest of the sections, we will follow the above prescriptions to calculate the e^, 
anti-proton, 7 and neutrino CRs in the U{l)B^3Li models. 



III. THE LEPTOCENTRIC U{1)b~3L, DARK MATTER MODEL 



Data from PAMELA show that there are excesses in cosmic energy spectrum, but not 



in anti-proton spectrum. This prompts several authors to propose leptophilic DM models 
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131]. The simplest leptophilic model is the gauged U{1)L,-Lj DM model [13|]. In the minimal 
SM, besides the SU (3)c x SU (2)^ x f/(l)y gauge groups, the family lepton number difference 
Li — Lj can be gauge without anomaly. The resulting Z' coupling to SM particles are 
leptophilic. Among the three possible models, two of the models L^. — and — can 
explain the PAMELA and ATIC data well, while the - model can fit PAMELA, 
FERMI and HESS data well and disfavor the other two models. This model leaving very 
little signal in hadronic mode. The anti-proton excess is suppressed at a negligible level for 
all ranges of anti-proton energies. 

With right handed neutrinos, one can gauge the B — L global symmetry without gauge 
anomalies. This model has also been studied in the context of the recent excesses. 
However, in this type of models, the Z' is non-leptophilic leading to sizable anti-proton 
excess in cosmic ray in contradiction with PAMELA data. 

To explain the PAMELA data, one does not need to have a pure leptophilic model. A 
modified leptophilic model with sufficiently large lepton fraction with non-negligible hadron 
fraction in DM annihilation can be a viable model. Because the not completely negligible 
hadronic fraction, there may be some interesting testable consequences. The leptocentric 



U{l)B~-iLi model [15| is an interesting example of this type. 

The t/(l)B_3Li model can be classified into three different models with Lj being: a) Lg, 
b) L^, and c) L^. Assuming the DM field ^/^ is a Dirac fermion which couples to Z' with a 
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charge a, the Z' couphng to SM fermions and DM are given by 



-'int " ^ 



3 

j=u,d,c,s,t 



(5) 



where /j and u^. are one of the lepton generations, z/r. is the right handed neutrino which can 
be hght (Dirac type of neutrino) or heavy (Majorana type of neutrino). To give a mass to 
Z', the simplest way is to introduce a SM singlet Higgs boson S with non-trivial U{l)B~3Li 
charge b. After the singlet develops a non-zero vacuum expectation value < S >= Vg/V^, 
the Z' mass mz> is given by m|, = iP'g'v'l. 

All of the three models a), b) and c) have suppressed hadronic DM annihilation fractions. 
This can be easily understood by noticing that the cross section of DM annihilate to quark 
final products a is proportional to 2 x 3 x 3 x (g' /3)'^ (the first factor 2 takes care of up and 
down sector, the first 3 is the number of quark generation, the second 3 is the number of 
color and last factor g' /3 is the Z' coupling to quark), but to charged lepton is proportional 
to {Sg'y. Therefore the ratio of final products with quarks to final products with charged 
leptons is given by 

Here we have assumed that the phase space with top quarks in the final states are approx- 
imately the same as other lighter quarks and leptons. This approximate is justified as the 
DM should have a mass of order TeV to cover the e"^ excess energy range. 

It is clear that this model has suppressed hadronic DM annihilation fraction, but non- 
zero at tree level. With the above suppressed hadronic fraction, all the above three models 
are allowed by the PAMELA data. But it is interesting to note, as will be shown later, 
that because the non-zero tree level cross section, just beyond the reach of the current 
PAMELA data these models predict noticeable anti-proton excesses. Therefore with higher 
experimental energy reach, these models will be tested. To further distinguish different 
models, one needs to consider more experimental constraints. 

Without carrying out detailed numerical analysis one can make some judgments on fa- 
vored models. Considerations from the e"^ spectrum shape can help to select preferred 
models. To this end we note that model a) would predict a sharp falling in energy 
around the DM mass because the DM annihilate directly into a pair of electron and po- 
sition through the Z' mediation. This model seems to be favored by the ATIC data, but 
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the more precise data from FERMI-LAT and HESS do not confirm the sharp falhng feature 
and therefore disfavored this modeL For models b) and c) the electron and position are 
secondary final products from fi~^fi~ (model b)) and r"*~r~ (model c)) after DM annihilation. 
The spectrum is much more softer and are favored by FERMI-LAT and HESS data, but 
can not fit the ATIC data well. 

Cosmic 7 ray can further distinguish models a), b) and c). Model c having r+r~ as the 
primary DM annihilation products has a large fraction in hadronic final state. There are a 
lot of vr*^ which can decay into 7 pairs which leads to visible 7 ray excesses at observable 
level. The present data from HESS disfavors this model. On the other hand, for models a) 
and b), 7 ray dominantly comes from final state radiation. We will show that 7 ray produced 
in this model is just below the current bound. 



IV. THE COSMIC RAY 



In the U{l)B~3Li niodel, the relic density of the dark matter is controlled by annihilation 
of tptp Z'* qjQj + lili + z/jZ/j, while the excess is mainly determined by tptp — 

Z'* lili. The interaction rate {av)i.i., with the k mass neglected is given by 

_ 9 a'g'^ml 

vr (s - m|,)2 + r|,m|, ' 

where v is the relative velocity of the two annihilating dark matter and s is the total DM 
pair energy squared in the center of mass frame. Tz' is the decay width of the Z' boson. If 
the Z' mass is below the ipip threshold which we will assume, the dominant decay modes of 
Z' are Z' Ylj QjQj + + i^ii^i, and Tz' is given by, neglecting final particle masses 

The annihilate rate {(Tv)reiic controlling the relic density is given by 31(o"f );.;- /IS. 

In the above, we have assumed that there are only left-handed light neutrinos. If there are 
light right-handed neutrinos to pair up with left-handed neutrinos to form Dirac neutrinos, 
the factor 31 in these equations should be changed to 40. 

Since the relic density of DM is determined by the annihilation rate {crv)reiic, the model 
parameters are thus constrained. The same parameters will also determine the annihilation 
rate producing the excesses observed today, which requires a much larger annihilation 
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rate. A boost factor in the range 100 to 1000 is necessary. We find that Breit-Wigner 
resonance enhancement mechanism works very well in our models if the Z' boson mass is 
about two times of the dark matter mass [ll|,[l3|,Q]. 

The boost factor in this case comes from the fact that since the Z' mass m^/ is close to 
two times of the dark matter mass m^, the annihilation rate is close to the resonant point 
and is very sensitive to the thermal kinetic energy of dark matter. To see this let us rewrite 
the annihilation rate into a pair of charged leptons as 

9c?g"^ 1 
^ \^Tim% (5 + t;V4)2 + 72 ' 

where we have used the non-relativistic limit of s = 4m^ + m^f^, with b and 7 defined as 



2 -4m2(l-5), and72 = r|,(l-5)/4m2. 



For thermal dark matter, the velocity is proportional to the thermal energy of dark 
matter. It is clear that for small enough b and 7, the annihilation rate is very sensitive to 
the thermal energy and therefore the thermal temperature T. At lower dark matter thermal 
energies, the annihilation rate is enhanced compared with that at higher temperature. This 
results in a very different picture of dark matter annihilation than the case for the usual 
non-resonant annihilation where the annihilation rate is not sensitive to dark matter thermal 
energies. The annihilation process does not freeze out even after the usual "freeze out" time 
in the non-resonant annihilation case due to the enhanced annihilation rate at lower energies. 
To produce the observed dark matter relic density, the annihilation rate at zero temperature 
is required to be larger than the usual one, and therefor a boost factor. With appropriate 8 
and 7, a large enough boost factor can be produced. For our numerical analysis we follow 
the procedures in Ref.[13|. 

As already mentioned in the previous section that if li = e, the is hard, but for k = fi 
or r, the have to come from fi and r decays leading to softer spectra. 

In the Figs. [H and [2] we show our results from fitting the PAMELA positron ratio and 
the ATIC or FERMI-LAT electron spectrum for the U{1)b-3l,, U{1)b-3l, and U{1)b-3l^ 
models taking the diffusion height Z^ 4kpc and 2 kpc respectively. 

There are several DM model parameters in the fitting, the DM mass mo, the gauge 
coupling constant ag'^. Since the Breit-Wigner enhancement for the boost factor forces the 
Z' mass rriz' to be 2mD, it is not an independent parameter in the models we are considering. 
We find that the DM mass for Lj = e, = /i and Lj = r are fixed to be in the ranges of 
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FIG. 1: The model prediction for cosmic ray energy spectrum. The upper panels show the 
positron ratio to fit the PAMELA data with t/(l)_B-3Li models. The lower panels show the electron 
spectrum to fit the ATIC or Fermi data. The height of the Galactic diffusion region is taken as 
4kpc. 
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FIG. 2: The same as Fig. [T] except the hight of the Galactic diffusion region taken as 2kpc. 



several hundred GeV to several TeV in order to cover the observed ranges of excesses. 
The corresponding values for a^f'^ are limited to be of order 10~^. The masses and the 
coupling ag'"^ used in the figures for three different models are 1 TeV, 1.5 TeV, 3 TeV, and 
1.7 X 10~^, 2.5 X 10"^ and 5.3 x 10~^, respectively. 
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To see how the astrophysics parameters affect the results, we show results for two different 
values for the hight of the diffusion region. In Fig. [T]we have taken the hight of the Galactic 
diffusion region as 4 kpc, while at Fig. [2] we take the hight 2 kpc. We have adjust the DM 
mass, the annihilation cross section to fit the data, and also the propagation parameters to 
account for all the cosmic ray data. For /j = e, that is for the U{l)B-3Le model, the 
spectra from DM annihilation are hard. With appropriate DM mass, from the figures we 
can see that this model can fit the PAMELA and ATIC data very well. For the other two 
models h = fi (the f/(l)B-3L^ model) or U = r (the ?7(1)b_3l^ model), the e"^ spectra are 
from secondary decays of /i or r and therefore are softer. It can be seen from the figures 
that these two models can fit the PAMELA and Fermi data. The U{1)b~3Lt model has an 
even softer spectra than those of the U{l)B-3Lf, model because r has large branching 
ratios into final states with no e in the final state. 



V. THE ANTI-PROTON COSMIC RAY 



The anti-proton from DM annihilation in the U{l)B-3Li comes from DM annihilation into 
quark pairs whose annihilation rate is given by 



qq 



Sn (s — m%,Y + T'i,mi, 



(10) 






FIG. 3: The ratio of anti-proton to proton as function of energy predicted in the B-3e, and 
B-3r models respectively from the left to the right panels with Zf^ = 4kpc. The cross section of 
DM annihilation is normalized to fit the ATIC or Fermi electron spectrum. 

The quark and anti-quark produced from DM annihilation are subsequently fragmentate 
into hadrons and may cause excesses in cosmic proton and anti-proton spectra. In our 
analysis the anti-proton spectrum from quark pair hadornization is treated with the PYTHIA 
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FIG. 4: The same as Fig. [3] but for Zh = 2kpc. 



package 26| . We find that different quark flavors have approximately the same proton and 
anti-proton spectra in the final products. Our results are shown in Fig. [3] and Fig. H] with 
the diffuse height Zh = 4kpc and 2kpc respectively. 

From the figures we see that at the high energy end of the anti-proton spectrum, there is 
already tension between the theoretical prediction and the current PAMELA data. For the 
diffusion hight Zh = 4kpc in Fig. [3] we notice that the antiproton fraction is obviously too 
much. However, due to large uncertainties in the propagation models, there are rooms to 
adjust parameters to fit the data. When we adopt the diffusion height Z^ = 2kpc, as shown 
in Fig. m the anti-proton flux is suppressed and the ratio is consistent with the PAMELA 
data at about 2a level. It should be noted that for both Z^. = 4kpc and Z^ = 2kpc we 
have adjusted the other propagation parameters, such as the diffusion coefficient D{E), so 
that all the cosmic ray data are well reproduced in both models. The DM annihilation cross 
section are normalized so that the ATIC or Fermi electron spectrum is well fitted, as shown 
in Figs. [Hand [21 

We also notice from Figs. [3] and [4] that the contribution of anti-protons from DM annihila- 
tion is most important at several hundred GeV. Therefore future data at higher energies will 
give much more stringent constraints on the DM contribution to anti-protons. We expect 
that the near future PAMELA and AMS02 data will play a crucial role in testing different 



DM annihilation models. 



VI. THE 7-RAY RADIATION FROM THE GALACTIC CENTER 



Since we require a very large DM annihilation cross section to account for the observa- 
tional CR excesses, this will lead to a strong 7-ray radiation from the final lepton states. 
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i.e. via ipip Z' ^ lihl- The cross section for high energy 7 of energy E^, with » ml^, 
is given by [2?! 



I") 
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(11) 

be 



27|. 



where s = 4m^ and s' = 4:m^{m^ — E^). For the r final states more 7-rays wil 
produced from the r hadronic decays, which will be calculated using PYTHIA package 
We will calculate the 7-ray emission from the Galactic Center (GC) in the models under 
consideration and check if they are consistent with the HESS observation at the GC 28 1. 

As can be seen from Eq. ([3]) that the 7-ray produced by DM annihilation depends on 
the square of the DM profile p and therefore the result will be sensitive to the form of p. In 
the following calculation we will take two popular dark matter density profiles, the NEW 



profile 



291 ] and the Einasto profile 30|], to predict the 7-ray fiux for comparisons. The two 



density profiles take the forms as 

pNFw{r) 

pEina{r) 



(r/rs) (1 + t/ts 
poexp ■ 



>2 ' 



2 

a 



(12) 



where we take the scale radius = 25kpc, a = 0.23 for the Einasto profile and normalize 
the local density as = 0.3GeV/cm^. Given the density profile, the 7-ray fiux along a 
specific direction can be calculated as given in Eqs. and (jl]). 

We calculate the diffuse emission at the GC with the galactic latitude b and longitude / 
in the region \b\ < 0.3°, |/| < 0.8° and compare the result with the HESS observation 28|. 
To compare with spectra reported by HESS we have followed the HESS collaboration to 
subtract the background from 0.8° < |6| < 1.5° [28 



In Fig. [5] we show the results together with the data taken from the HESS observation 
]. For the r final states we notice the 7-ray fiux from GC is greatly 



on the same region 



higher than the observed one. Therefore the B — 3t model is disfavored considering its 7-ray 
emission at the GC. Further, 7-rays also exceed the HESS result if taking the NEW profile 
no matter what the final states are taken. However, the recent high precision simulation 
show that the density profile at the halo center does not tend to any simple form of power 



law [31|- Instead the Einasto profile gives best fit to the most recent simulations |31l . I32 |. 
We also notice that if DM annihilates into electrons or muons the 7-ray fiux is consistent 
with HESS when taking the Einasto profile. 
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FIG. 5: The 7-ray emission from the GC reg ion in the B — 3Lj models. The data are taken from 



28(1 . We have taken NFW and Einasto profiles to predict 



the HESS observation at the same region 
the 7-ray fluxes. 

It should be noted that for the two propagation models with = 4kpc and 2kpc the 
7-ray flux from the GC is almost the same, since the HESS observation focus on a very small 
region around the GC. 



VII. THE NEUTRINO COSMIC RAYS 

In the U{l)B-3Li models, there are large fractions that DM annihilate into neutrino and 
anti- neutrino UiUi pairs directly. The annihilation rate is given by 

9 



(13) 



vr (s — m|,)2 + r|,m|, 

where A = 1 if the right- and left-handed neutrinos pair up to have Dirac mass, and A = 1/2 
if the right-handed neutrinos are heavy. This is a two body annihilation process and the 
neutrino energy spectrum produced will be line spectrum type. Therefore the line spectrum 
neutrinos, if detected, are distinctive signals of DM annihilation, since no astrophysical 
processes can produce such line spectrum. 

There are also neutrinos from secondary decays of fi ^ i^^iCUe for the ?7(1)b-3Lp and 
T z/T-/iz/^ + i^^ei^e for the U{1)b~3l^ model, r can also produce softer neutrinos by hadronic 
decays. However, these continuous spectrum neutrinos are hard to be discriminated from 
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the atmospheric neutrino background. Further, taking these soft neutrinos into account we 
have to lower the threshold of the neutrino detector. The amount of atmospheric neutrinos 
increase rapidly by lowering energy since they have a very soft spectrum 33 1. 

There have been some studies of detectability of neutrinos from the GC by DM annihi- 
lation 3J]. As shown in the last section the B — 3t model has been excluded by the 7-ray 
emission. Therefore we calculate the neutrino fluxes in the cases of -B — 3e and B — 3fi with 
the Einasto profile which are consistent with the HESS 7-ray observation. We would like to 
check if the future IceCube observation can detect the neutrinos predicted by our models at 
the GC. 

We will consider the muon neutrino detection rate at IceCube. The flavor eigenstates 
produced at the source will be modified at the detection point due to neutrino oscillations. 
The relation between different flavors from the GC to the Earth is given by 351] 



H Earth 



ijWj)GC , 



(14) 



where Pij is the probability of converting the j-th flavor type of neutrino at the source to 
the i-th flavor neutrino at the detection point. 

Assuming tri-bimaximal neutrino mixing pattern, the Pij matrix after neutrinos traveling 
a galactic distance is given by 

'^°-4cu 2cu 2uj ^ 

2u i-uj 4:-uj , (15) 
2cu A-cu A-cu 



1 



where ui = sin^(2^i2) ~ 0.87. Since the tri-bimaximal mixing is a good approximation to 
present data, we will use the above formula for numerical estimate. We find that the muon 
neutrino fluxes arriving at the Earth are 0.22 and 0.4 times of the initial and i'^ fluxes at 
the sources, respectively. 

We consider the detectability of the line neutrino spectrum of our models by the planned 
additional DeepCore within the IceCube detector. We take a conservative angular resolution 
as large as vr Sr. Considering the contained events within the DeepCore it is not difficult to 
calculate the z/^ events for 4 years observation with the DeepCore volume 0.04km^. The total 
contained muon events are about 606 for z/^ initial flavor coming from the vr solid angle cen- 
tered at the GC when assuming the DM annihilation cross section is large enough to account 
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for the cosmic e"^ anomalies and the Einasto density profile. Taking the energy resolution of 
the detector ailog-^^Q^E)) ~ 0.4 [36| (corresponding to the energy bin GOOGeV ~ 3.77TeV) 
and taking the atmospheric neutrino background of 3^ we get about 2350 contained back- 
ground events at DeepCore/IceCube for 4 years observation from the GC direction for n Sr 
solid angle. About 68% signal events will fall within the la energy range. Then we find that 
the DeepCore/IceCube will be able to detect the line spectrum neutrinos at the confidence 
level 8.5cr signal for 4 years observation with the B — 3fi model, if the DM annihilation rate 
is fixed by the excess data. For the B — 3e model (with the u^. initial flavor) the DM 
mass is ITeV. Therefore the energy bin is from AOOGeV to 2.5TeV. We then get the moun 
contained events for 4 years observation are about 243. The background events are much 
larger with about 4000 events since the energy bin is lower now. In this case we can only 
have a 2.6a signal for the z/g initial state which will oscillate into i^^. 

Here we have taken a conservative angular resolution. Considering that the DeepCore 
may have a better angular resolution we may get signals with much higher significance. 
Therefore the neutrino events are very promising signals to test the DM annihilation model 
which explain the cosmic excesses. If the scenario is correct we expect to detect at 
DeepCore/IceCube of the line spectrum neutrinos predicted in the B — 3Lj models. 



VIII. CONCLUSIONS 

In this work we studied cosmic e^, p, 7 and neutrino rays in the leptocentric ?7(1)b_3l. 
dark matter models. In these models, DM annihilation into SM particles is mediated by 
the Z' gauge boson of U {1)3-31,- The couplings of Z' to leptons are larger than those to 
quarks leading to larger annihilation of DM to leptons than hadrons. This naturally explains 
the excesses from PAMELA, ATIC or FERMI observational data, and at the same time 
these models can suppress the anti-proton flux efficiently. We find that the U{1) b-sLe model 
can fit PAMELA and ATIC, while U{1)b-3^ and U{l)B-3r can fit PAMELA and FERMI, 
respectively. Future data can be used to distinguish different models. We showed that 
by slightly adjusting the propagation parameters these models can predict anti-proton flux 
consistent with the PAMELA data. But these models are different than pure leptophilic 
models which predict negligible anti-proton cosmic ray in all energy ranges, the leptcentric 
U{l)B-3Li models predict anti-proton flux beyond the background at higher energies which 
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can be tested by the near future data from PAMELA or AMS02. 

The U{l)B-3Li niodels have predictions for cosmic 7-rays from DM annihilation. We 
investigated the 7-ray emission from the GC in these models. We found that the U{1)b~3t 
model is excluded by the HESS observation of the GC region since this model predicts too 
much 7-rays. Adopting the Einasto DM density profile the other two models predict 7-rays 
consistent with HESS data. 

We also investigated the detectability of neutrinos from the GC by the U{l)B-3e and 
U{1)b-~3ij, models at DeepCore/IceCube. In these models the fraction of the DM annihi- 
lation into neutrinos are large and produce line spectrum neutrinos detectable by Deep- 
Core/IceCube. We found that the DeepCore/IceCube can have a 8. 5(2. 6)0" signal of the 
line neutrino spectrum for 4 years data taking. This will be a distinctive signal of DM 
annihilation in our leptocentric U{l)B-3e and U{l)B^3fi models. 

Finally we would like to have some comments on the detectability of leptocentric DM 
model at the LHC and direct DM detection experiments. Since the leptocentric DM models 
discussed here are very different than the pure leptophilic DM models, where the Z' only has 
interactions with leptons at the tree level, the Z' in our cases has interaction with quarks 
and may lead to detectable effects at the LHC and also direct DM search on the Earth. 
However because we rely on Breit-Wigner resonant enhancement effect to produce the large 
boos factor, the Z' mass is constrained to be 2 times of the DM mass and is constrained to 
be in the TeV range and the couplings ag'^ to be in the range of 1 x 10~^ to 5 x 10"^, the 
cross section for producing Z' and virtual effects due to Z' for the LHC and also for direct 
detection are small, it is not possible to have observable effects at the LHC and near future 
direct DM detection experiments. To have a lower Z' mass and therefore to have observable 
effects at the LHC and direct DM detection experiments, a different mechanism for the 
boost factor has to be in effective to relax the relation that the Z' mass is about 2 times of 
the DM mass. To this end we note that if there is a long range interaction between DM by 
exchanging a light new particle, it is possible to have the Sommerfeld effect in operation, 
e.g. m Ref. HQ- 

This offers another interesting possibility to explain the cosmic data. 
We will give details for this possibility elsewhere. 
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